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Abstract 

We calculate the quantum phase diagram of the XXZ chain with nearest- 
neighbor (NN) Ji and next-NN exchange J2 with anisotropies Ai and A2 
respectively. In particular we consider the case Ai = — A2 to interpolate 
between the XX chain (Aj = 0) and the isotropic model with ferromagnetic 
J2. For Ai < —1, a ferromagnetic and two antiferromagnetic phases exist. For 
|Aj| < 1, the boundary between the dimer and spin fluid phases is determined 
by the method of crossing of excitation spectra. For large J2/ Ji, this method 
seems to indicate the existence of a second spin fluid critical phase. However, 
an analysis of the spin stiffness and magnetic susceptibility for Ai = A2 = 1 
suggest that a small gap is present. 
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I. INTRODUCTION 



In recent years, several interesting quasi one dimensional magnetic systems have been 
studied experimentally [jl]-ffl. Among them some compounds containing CuO chains with 
edge sharing CUO4 plaquettes are expected to be described by the XXZ model with an 
important next-NN exchange J 2 in comparison with the NN exchange Ji, and the latter 
can also be negative if the Cu-O-Cu angle is near 90° Q. In addition, there has been 
a renewed interest in this XXZ model (which is equivalent to a zig-zag ladder) recently 
P-IB|. In particular, for the XX chain (Ax = A 2 = 0), field theoretical methods predicted 
a critical (gapless) phase with incommensurate correlations for J 2 >> J\ > 0. This seems 
confirmed by recent numerical work JTof- The effects of magnetic field and magnetization 
jumps (metamagnetism) have also been studied recently [pl]-14|. Metamagnetism is more 
likely to occur in nature for ferromagnetic Ji, because smaller anisotropies are required 
|T3| . Interestingly, for ferromagnetic small J\ a very small gap with an astronomically large 



associated correlation length has been found by field theory calculations and seems consistent 
with density matrix renormalization group (DMRG) calculations ||16|| . 

The phase diagram of the model for Ai = A2 > has been accurately determined 
by Nomura and Okamoto using the method of level crossings of the excitation spectra, 
supported by results of conformal field-theory and renormalization group |T7| . This method 
has been also used successfully to determine phase diagrams of electronic models |1| , and is 
related in some cases to jumps in topological numbers determined by Berry phases fll9l , |20" 



Using the method of level crossings, the phase diagram has been extended recently to the 



region Ai = A 2 < fL5 |. However, to assume that A 2 < implies a different sign of the z 
component of J 2 with respect to those of the x, y plane. This is an unrealistic anisotropy (this 
inconvenient is not present for J\ because the sign of its x, y components can be changed 
rotating every second spin in the chain by it around the z axis). In addition, the above 
mentioned systems containing CuO chains with edge-sharing plaquettes and Cu-O-Cu angle 
near 90°, are expected to lie near the isotropic case A 1 = — 1, A 2 = 1. 



2 



In this work we calculate the phase diagram of the system for A2 = |Ai| and —1.2 < 
Ai < 1.2 using the method of level crossings. This includes the interpolations between the 
XX chain and the isotropic cases with ferro- or antiferromagnetic J\. 



II. MODEL AND METHOD OF LEVEL CROSSINGS 

The model is: 

H = / ][Ji(SfSf +1 + SfS^ +1 + A 1 S-S- +1 ) 

i 

+J2(SfSf +2 + S i S i+2 + A 2 SiSf +2 )], (1) 

where S? is the (3 component of the spin- 1/2 operator at site i. 

The basic idea of the method of level crossings is that the properties of a critical phase 
at low energies and large distances are ultimately determined by a scale invariant fixed 
point [|17|. Then, using conformal field theory one can relate the excitation energy which 
corresponds to some operator A4 at site % (for example a spin flip S~), to the dependence 
of the correlation functions of this operator with distance: 

E A (L)-E g (L) = 2 -^, (A +d A,)^J^. (2) 

Here L is the length of the system, v the spin- wave velocity, E g (L) the ground state energy, 
and xa the critical dimension for the excitation A. For example if J 2 = 0, the XXZ model 
was exactly solved and its properties are well known |2l|j22| . The long distance dependence 
of the main correlation functions is 

(S&Sr) - (-V d r- 1/K , (S? +d S*) - {-l) d r- K (3) 

with 

K = ~ / a \ (4) 

7r — arccos(Ai) 

if I Ai| < 1. For I Ai| < 1, the minimum excitation corresponds to a spin flip with wave vector 
71. The ground state has total spin projection S z = ±1, its wave vector differs from that of 



the ground state in tt, and its parity under inversion P is opposite to that of the ground state 
fTTfl . The matrix element of S^~^ between the ground state and the first excited doublet 
is different from zero. Then, the correlation functions (S^~ +d S^ are the dominant ones at 
large distances. 

At Ai = 1, the above mentioned doublet crosses with the excited state with S z = 0, 
(with wave vector differing from the ground state in tt, and P and parity under spin reversal 
T opposite to that of the ground state). Then, for Ai > 1, (S z +d S z ) are the dominant 
correlations at large distances. At Ai = 1, a gap opens and conformal invariance ceases to 
be valid, but for a sufficiently small gap so that correlation length is much larger than the 
size of the system, Eqs. (0) are still expected to hold. 

In summary inside a critical phase or near its boundaries, the quantum numbers of the 
first excited state determine the dominant correlations at large distances, and therefore the 
crossings of the excited states determine the boundaries. 

For small J2/J1, the use of renormalization group allows to determine the size dependence 
of the crossings and to extrapolate them very accurately to the thermodynamic limit jTTJ] . 



III. RESULTS 

The phase diagram we obtain in the plane a = J^j J\ as a function of A 1; with A 2 = |Ai| 
is represented in Fig. 1. The transition between dimer and spin fluid I phases for A x > 



has been obtained before by Nomura and Okamoto |T7|]. The spin fluid I phase is the well 
known ground state of the model for a = 0, |Ai| < 1. The dimer or spin Peierls phase has a 
spin gap and broken translational symmetry (the unit cell is doubled) in the thermodynamic 
limit. The first excited state in finite systems has the same quantum numbers P and T as the 



ground state, but its wave vector differs by n |L7|]. The new results for the spin fluid I-dimer 



transition are displayed in Fig. 2. The transition points (open squares) were accurately 
determined extrapolating the crossing of the first excited state using data for chain lengths 
L = 12, 16 and 20. As explained earlier |T7| the crossing point has a 1/L 2 size dependence. 
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We have verified this also for negative Ai (see Fig. 3). However, in the region denoted by 
? in Fig. 3 (constructed with data for L = 16), the first excited state has opposite P and 
T as the ground state. This is unexpected for the dimer phase. This might be a finite-size 
effect related with the high degeneracy at the point Ai = — A2 = —1, a = 0.25 fll3 , p5| . The 
other possibility is that a novel phase exists there, but this seems unlikely. 

After the transition to the dimer phase, if a is increased further, the crossing of excitation 
spectra indicates another transition to a phase in which the first excited state has S z = ±1 
and wave vector k c = it/ 2. We denote this phase as "spin fluid II " (between quotation 
marks), because, as we discuss below, its nature is not established for all values of Aj. For 
Aj = 0, from field theory results one might expect a spin fluid phase for large a ||. One also 
expects that k c is actually incommensurate and slightly different from tt/2 @J^J^,§]. However, 
we are not able to detect this difference using periodic boundary conditions. While the use of 
twisted boundary conditions has given interesting results concerning this incommensurability 
||, to implement them in the level crossing method is time consuming and present technical 
complications, which we avoid in this work. 

As shown in Fig. 4, the level crossing for the dimer-spin fluid II transition does not follow 
a 1/L 2 dependence, and the extrapolation to the thermodynamic limit is not so accurate as 
for the dimer-spin fluid I transition. We have used a quadratic extrapolation in 1/L 2 , using 
the data for L — 8, 12, 16 and 20. Field-theoretical results suggest that there is a small gap 
and therefore no spin fluid phase for Ai = A 2 = 1, and J 2 >> J\ HJ7|Jl^|. For A x = —1, 



A 2 = 1, a tiny gap is predicted with an associated astronomically large length scale |]T6 
A direct extrapolation of the gap from our finite-size results has an error larger or of the 
order of the gap itself and therefore, it is not able to establish if the gap is open or not. 
However, in principle we can check if a given phase is a spin fluid phase by other methods. 
For example, the energy per site should vary with system size as: 

e(L) = e(oo) - -^v s (5) 

where v s is the spin velocity (or the sum of spin velocities if there were several types of 



low-energy excitations). We have checked that Eq. ([|) is not satisfied inside the dimer 
phase. A fit of the energies for L = 12, 16, 20 and 24 with Eq. (|5|) give an error of v s which 
is of the order of v s itself. Instead, inside both "spin fluid" phases the fit is good. The error 
in v s is of order of 0.1% inside spin fluid I, and of order of 1% inside spin fluid II. Within 
spin fluid I, we also find that several relations derived from conformal field theory hold, and 
are consistent within them within ^ 5%. The spin velocity can be calculated as: 

E(2n/L) - E(0) 
Vl = 2^1 ' (6) 

where E(q) is the lowest energy in the S z = sector. The extrapolation of v\ to the ther- 
modynamic limit using a quadratic polynomial in 1/L is only slightly less than v s obtained 
fitting Eq. (||). Also from the numerical calculation of the susceptibility x an d spin stiffness 



XS B L[E(1) + E(-1)-2E{0)] (7) 
^ Ld 2 E 

where here E(S Z ) is the lowest energy for total spin S z , and is a flux opposite for spin up 



and down, B we can obtain K in two independent ways from the relations p2 



The correlation exponents obtained from these two equations and Eqs. @, (0), @ with 
Xi — Xai D\ — D s also agree within 2% to 10% in the spin fluid phase I. The agreement 
improves with increasing Ai. 

A simple analysis based on Eqs. ([?]), (§), @ is not possible in the "spin fluid II" phase. 
Eqs. @ and @ suggest that the low energy properties of this phase are given by two free 
bosonic theories (this is of course true if J\ — 0): v s — V\ +V2, with v% slightly larger than v\ 
(see Fig. 5(a)). This is in principle reasonable, since if there are two branches of low-energy 
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excitations, Eq. @ should give the smallest velocity. Also one expects that Eq. (§) gives 
the sum of both spin stiffness: 

D s = KlVl + K ^ (10) 
An 

The case of the susceptibility is more delicate. From thermodynamics, the total susceptibility 
should be the sum of those of both branches. However, our numerical study suggests that 
when one flips only the spin, it goes to the branch of lowest velocity. 

The results of field theory in the case J 2 3> J\ |f7|J§ (and also those of the Hubbard ladder 
for small interchain hopping |24|,25|) suggest that the effective low-energy theory splits in two 



sectors, symmetric tp + and antisymmetric </?_, and in general, at least the latter is gapped. 
For Aj = A 2 = 0, ip + is massless and y?_ has a gap which depends on J2/J1 as a power 
law |||9|]. In the isotropic case A x = A 2 = 1, both sectors have an exponentially small 
gap 0. These field theory results are inconsistent with the hypothesis of two branches of 
massless excitations in the "spin fluid II" phase. To check this hypothesis we investigated 
the consistency of Eqs. (|5|) with v s = v\ + v%, (^j), ( |TD| ) and an expression for the total 
susceptibility, for some values of K\ and K 2 in the above mentioned cases (A x = A 2 = 0, 
and Ai = A 2 = 1). We do not find reasonable results. For example, eliminating Ki from 
these equations, sometimes one negative result arises. The most plausible explanation is 
that at least one of the modes is gapped, but the system is not large enough to detect it 
in the size dependence of the energy (Eq. (|5]) is only valid for sufficiently large L). In fact, 
some of our fittings suggest that the magnitude of the slope of e vs 1/L 2 is decreasing with 
system size. 

To end this section, we discuss the observed phases for 1.2 > |Aj| > 1. For Ai < —1, 
the ground state changes from a fully polarized ferromagnet (S z = L/2) to S z = at 
the solid squares of Fig. 1. Actually, there is a small region of intermediate S z between 
both phases, but it decreases with system size an seems to disappear in the thermodynamic 
limit. In the phase with S z = 0, the first excited state has wave vector tt/2. This fact, 



a study of correlation functions |26|] and a trivial analysis of energies in the classical limit 



— Ai = A 2 — > oo, suggests that the phase has long-range AF order of the type TTIITT ••• 
(denoted as AFII). The same features denoting the presence of the AFII phase, are also 
found for A; > 1 and large a. We expect that for |Aj| > 1, a spin wave approximation can 
describe qualitatively the essential physics. 



IV. DISCUSSION 

Using the method of level crossings, we have calculated the phase diagram of the XXZ 
model with next-NN exchange (Eq. (j|D) for A 2 = |Ai|. For the dimer-spin fluid I transition 



our results are quite robust and extend previous results |ll7||l5| to the case 1 < Ai = — A 2 < 
0. In this region, the transition takes place near a ^ 0.3. For larger values of a the method 
predicts a transition from the dimer phase to a second spin fluid phase. The transition is at 



a c = 1.42 for A x = A 2 = 0, in qualitative agreement with a c ^ 1.26 obtained in Ref. |T0 
For these values of Aj, field theoretical results predict the existence of this phase ||. 

For Ai = —1, a tiny gap with an astronomically large length scale is predicted in Ref. 
|T6|j . Such a length is of course much larger than the system sizes we use, and our finite size 



results are consistent with those of Ref. [[L6[| . 

For Ai = A 2 = 1, we have tried to interpret the thermodynamic properties of the "spin 
fluid II" phase as a sum of two independent fluids, but the results are not consistent. We 
conclude then that at this point, the extrapolation of the level crossings is not reliable and 
in agreement with other calculations ||6|,|7|Jl4]|, the gap persists for large J 2 . 

From the parameters estimated for several compounds near the isotropic limit Ai = 
-A 2 = -1 : La 6 Ca 8 Cu 24 4 i (a - 0.36), Li 2 Cu0 2 (a - .62) and Ca 2 Y 2 Cu 5 Oi (a - 2.2) || 
our results indicate that these compounds lie outside the usual spin fluid I phase. They are 
in the AF2 phase if Ai < —1. If Ai > — 1, LagCasCu^O^ is clearly in the dimer phase. One 
expects a very small, perhaps unobservable gap for Ca 2 Y 2 Cu50io and spin-spin correlations 
with wave vector near n/2. Li 2 Cu0 2 lies in between. The interchain interactions might 
affect this scenario. 
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FIGURE CAPTIONS 

Fig.l. Phase diagram of the model Eq. (P in the Ai, a = Jij J\ plane, keeping 
A 2 = | Ax | . F denotes the fully magnetized ferromagnetic phase, and AF2 a phase with long 
range order fill ••• • 

Fig. 2. Details of the boundary of the spin fluid I phase of Fig.l for negative Ai. 

Fig. 3. Size dependence of the crossing between dimer and spin fluid I excitations. 

Fig. 4. Size dependence of the crossing between dimer and "spin fluid IF excitations. 
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